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ABSTRACT 

MCG-5-23-16 is one of the first active galactic nuclei (AGNs) where relativistic reverberation in the iron K line 
originating in the vicinity of the supermassive black hole was found, based on a short XMM-Newton observation. In 
this work, we present the results from long X-ray observations using Suzaku , XMM-Newton, and NuSTAR designed 
to map the emission region using X-ray reverberation. A relativistic iron line is detected in the lag spectra on three 
different timescales, allowing the emission from different regions around the black hole to be separated. Using 
NuSTAR coverage of energies above 10 keV reveals a lag between these energies and the primary continuum, which 
is detected for the first time in an AGN. This lag is a result of the Compton reflection hump responding to changes 
in the primary source in a manner similar to the response of the relativistic iron K line. 
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1. INTRODUCTION 

The primary X-ray continuum emission in active galactic 
nuclei (AGNs) is due to Compton scattering of accretion 

disk photons in a hot corona (Haardt & Maraschi 1991). 
The relatively cold and dense surrounding material intercepts 

and reflects some of this emission (George & Fabian 1991), 
producing a characteristic reflection spectrum (Ross & Fabian 

1993; Garcia et al. 2013). Relativistic effects imprinted on the 
reflection spectrum provide a measure of the distance of the 

emitting region from the black hole, allowing its spin to be 

directly measured (Fabian et al. 1989; Laor 1991; Brenneman 

& Reynolds 2006; Miller et al. 2009; Reynolds 2012). The fast 
variability in the primary emission (e.g., Vaughan et al. 201 1) is 

also mirrored by the reflecting medium (through a reverberation 
transfer function), providing direct measures of physical scales 
and geometries (Matt & Perola 1992; Campana & Stella 1995; 
Reynolds et al. 1999; Wilkins & Fabian 2013). 

After a hint of detection in Ark 564 (McHardy et al. 2007), 
reverberation lags were significantly observed in several Seyfert 
galaxies, where the soft (0.5-1 keV) band, composed mostly of 
reflected iron L emission and other species, is observed to lag 
the harder (1-4 keV) band, which is dominated by the primary 
emission (Fabian et al. 2009; Zoghbi et al. 2010; Zoghbi & 
Fabian 2011; Emmanoulopoulos et al. 2011; de Marco et al. 
2011; Fabian et al. 2013; Kara et al. 2014; Cackett et al. 2013). 
The time delay and the timescale on which it is observed both 
appear to scale with black hole mass (De Marco et al. 2013). 

A major step forward was the discovery of iron K reverber- 
ation in NGC 4151 (Zoghbi et al. 2012). Whereas the earlier 
lags were between the primary continuum and a blend of rel- 


ativistically broadened emission lines (comprising the soft ex- 
cess), the lags in NGC 4151 were between the same primary 
continuum and the iron K line at ~6.4 keV. The significance 
of the iron K reverberation is that it is a single line in a clean 
part of the spectrum, allowing direct spectro-timing modeling 
that directly probes the emitting region (Cackett et al. 2014; 
Emmanoulopoulos et al. 2014). Iron K line reverberation has 
now been observed in several other objects (Zoghbi et al. 2013b; 
Kara et al. 2013b, 2013a; Marinucci et al. 2014), and again the 
delay and the timescale on which it is observed, both appear to 
scale with mass (Kara et al. 2013b). 

In this work, we study the reverberation in MCG-5-23-16 
in detail and present the first measurement of a lag in the 
Compton hump in an AGN. A new observation of the source with 
the Nuclear Spectroscopic Telescope Array ( NuSTAR ; Harrison 
et al. 2013) revealed a time delay above 10 keV that is consistent 
with reverberation in the innermost regions of the accretion disk. 

MCG-5-23-16 is a Seyfert 1.9 Galaxy (z = 0.0085) with 
a typical 2-10 keV flux of ~8 x 10 1 1 erg cnT 2 s _1 and a 
mass of 1-5 x 10 7 M e (Wandel & Mushotzky 1986; Ponti 
et al. 2012). Its spectrum resembles a classical Compton-thin 
Seyfert 2 galaxy with a column that does not affect the spectrum 
above 3 keV significantly. The spectrum below 1 keV contains a 
combination of emission from scattered continuum photons and 
distant photoionized gas. Above 2 keV, the spectrum shows both 
narrow (EW-60 eV) and broad (EW- 50-200 eV) iron K a 
lines along with a strong Compton hump above 10 keV (Weaver 
et al. 1997; Balestra et al. 2004; Mattson & Weaver 2004; Braito 
et al. 2006; Reeves et al. 2007; Zoghbi et al. 2013b). 

This work uses new observations of MCG-5-23-16 with 
Suzaku, XMM-Newton and NuSTAR with three specific aims: 
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Figure 1 . 2-10 keV light curves from all the new observations of MCG-5-23-16. Top: Suzaku XIS (combined front and back-illuminated) from two observations 
(ObsIDs: 708021010 and 708021020), shown in blue (left axis). The panel also show in orange (right axis) the simultaneous NuSTAR light curve (ObsID: 60001046002), 
where data from FPMA and FPMB have been combined. Bottom: XMM-Newton EPIC PN light curves from the full two orbits (ObsIDs: 0727960101 and 0727960201, 
not simultaneous with the top panel) taken three weeks later. 

(A color version of this figure is available in the online journal.) 


1 . To confirm the iron K reverberation lags reported in Zoghbi 
et al. (2013b) using a single XMM-Newton observation. 

2. To use the combined data sets to measure lag-energy 
spectra on three timescales, localizing the emission to three 
emission regions 

3. To observe time lags due to reverberation in the Compton 
reflection hump using NuSTAR. 

2. DATA REDUCTION AND ANALYSIS 

In order to study the detailed variability of MCG-5-23-16, it 
was targeted with a long Suzaku observation starting on June, 
01 2013 for a net exposure of 300 ks (ObsIDs: 708021010 and 
708021020). NuSTAR observed it simultaneously starting on 
June 03, 2013 for an exposure of 150 ks (ObsID: 60001046002). 
XMM-Newton (Jansen et al. 2001) also observed MCG-5-23-16 
for two full orbits three weeks later starting on June 24, 2013 
(ObsIDs: 0727960101 and 0727960201). We concentrate in this 
work on the timing results. Spectral modeling will be published 
separately. 

Suzaku and NuSTAR data were reduced using heasoft ver- 
sion 6.15 with the latest calibration files (CALDB files released 
December 23, 2013 and February 3, 2014 for NuSTAR and 
Suzaku, respectively). XMM-Newton data were reduced us- 
ing sas ( xmmsas_20131209_1901-13.5.0 ) and also using the 
latest ccf calibration files including the EPIC-pn (Striider et al. 
2001) long-term charge transfer inefficiency corrections. Suzaku 
light curves were extracted from different energies (Section 3) 
using xsefect to filter events from a source and background re- 
gions. Background counts scaled to the source region size were 
then subtracted from the source counts. NuSTAR light curves 
for source and background regions were corrected for livetime, 
psf, exposure, and vignetting using the task nulccorr. Abso- 
lute and relative corrections for the EPIC PN and MOS light 
curves were applied using epiclccorr. For the purpose of the 
analysis presented here, all light curves have time bins of 5 12 s. 
Other time bins do not change the results. 

Time lags are calculated using the method detailed in Zoghbi 
et al. (2013a). The method uses maximum likelihood to ob- 
tain the frequency-resolved time lags by directly fitting the light 
curves (Miller et al. 2010). The uncertainties in the lag measure- 
ments represent the 68% confidence intervals calculated as the 


values that change — 21og(£/£ max ) by 1 where C is the likeli- 
hood value and £ max is the optimum likelihood value (see details 
in Zoghbi et al. 2013a). This means the lag errors in the plots 
represent lcr uncertainties. The errors quoted in the text when 
doing model fitting represent the 90% confidence limits for one 
interesting parameter (i.e.. A/ 2 = 2.71). 

3. RESULTS 

3.1. Iron K Reverberation 

In this section we concentrate on the 2-10 keV energy band. 
Here we use the combined XMM-Newton, Suzaku, and NuSTAR 
(< 10 keV) data sets to study time lags in the iron K line. We first 
present a simple method for directly observing lags in the light 
curve. We then present a more general analysis using frequency- 
resolved time lags. 

3.1.1. Simple Considerations 

For a start, we consider a simple procedure to study the lag as 
a function of energy. The light curves from the new observations 
are plotted in Figure 1 . The top panel shows the simultaneous 
Suzaku XIS (longest blue) and NuSTAR light curves (orange), 
while the XMM-Newton EPIC-PN light curve is plotted in 
the lower panel. The light curves are characterized by high 
variability and strong flares. 

It is known from a previous XMM-Newton observation 
(Zoghbi et al. 2013b) that the peak of the iron line (~6.4 keV) is 
delayed with respect to lower and higher energies on timescales 
longer than ~ 1 0 ks. The simplest test one can do is look for 
lags within flares of a particular timescale at different energies. 
Although most flares contain multiple timescales (i.e., small 
scale flares on top of a longer timescale trend), the single flare 
at the end of the first XMM-Newton orbit (just below 10 5 s in 
Figure 1 (bottom)) is particularly interesting. It appears mainly 
as a ~10 ks flare. 

Figure 2 (left) shows a zoom in of this flare from eight energy 
bins between 2 and 10 keV in steps of 1 keV. The flare in each 
energy band was fitted with a simple Gaussian model to obtain 
a measure of the peak of the flare. The best fitted Gaussian 
centroids and their widths are plotted in Figure 2 (right). The 
best-fit Gaussian peaks were shifted vertically with the same 
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Figure 2. Left: a plot of the light curve section showing the strong flare at the 
end of the first XMM-Newton light curve plotted at different energies (energy 
increases down, and the x-axis has been shifted to zero). The lines on top of each 
flare represent the best fitted Gaussian model (see Section 3.1.1). Top right: the 
best fit values for the peak of each Gaussian model from the left panel (i.e., the 
lag of each flare compared to the 2-3 keV flare) plotted as a function of energy. 
There is a clear feature at the 6-7 keV band. Bottom right: the best fit values for 
the width of each Gaussian model from the left panel. The width of the flares 
also changes with energy, and tracks the shape of an iron line similar to the flare 
peak. 

(A color version of this figure is available in the online journal.) 


amount so that the first point has zero lag. It is clear from the 
change of the peak of the flare that the iron line lags the bands 
above and below it. The shape of the “broad iron line feature” is 
strikingly similar to that obtained from full Fourier analysis of 
the previous observation (Zoghbi et al. 2013b) despite the fact 
that it is extracted only from a single flare. Figure 2 (right) also 
shows the width of the flare as a function of energy. Interestingly, 
this too has a broad-line-like feature peaking at 6-7 keV. 

3.1.2. Frequency-resolved Time Lags 

The previous method does not take into account the stochastic 
nature of the light curve. Therefore, in order to study the lags in 
a systematic manner, using all the available data, we explore the 
frequency-resolved time lags. Although the XMM-Newton PN 
light curves are evenly sampled, those from Suzaku and NuSTAR 
are not due to Earth occultations. We use a maximum likelihood 
method to take this into account by directly fitting for the time 
lags at different temporal frequencies. The method divides the 
observed frequency range into a number of bins taking the lag 
value at each bin as an unknown parameter. A covariance matrix 
is constructed by an inverse Fourier transform and is used to 
calculate a likelihood function. Assuming uniform priors on the 
parameters and using a Bayesian formalism, lags at different 
frequencies and their probability distribution are obtained (see 
Zoghbi et al. 2013b). The end result when comparing two light 
curves is a lag-frequency plot. If more than one light curve 
is available (e.g., from different energies), each light curve is 
compared to a reference light curve, and the lag-energy plot can 
be constructed for the temporal frequencies of interest. 

Figure 3 shows the lag-frequency plot for the peak of the 
iron line (6-7 keV band) compared to both lower (2-3 keV; 
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Figure 3. Frequency-resolved time lags for MCG-5-23-16 comparing the peak 
of the iron K line to its wings from combined analysis (first row) and from 
individual instruments (rows 2-4, as labeled). The left column compares the 
6-7 keV band to the 2-3 keV band (blue circles). The right column (orange 
triangles) compares the 6-7 keV band to the 9-10 keV band (for XIS, we 
use 7-8 keV instead of 9-10 keV to obtain a better signal). We follow the 
usual convention, where positive means hard lags, i.e., hard band lagging soft 
band. The left column shows that the 6-7 keV band lag bands below it (positive 
lag). The right column shows that the 6-7 keV lags bands above it (negative 
lags). The different rows show that this conclusion is consistent between the 
different instruments used. 

(A color version of this figure is available in the online journal.) 


left column) and higher energies (right column). Light curves 
from the different detectors were fitted together by constructing 
a joint likelihood as well as separately. Positive and negative 
values indicated hard (i.e., hard bands lag soft bands) and soft 
lags, respectively. Two clear results can be seen in Figure 3: 

1 . The 6-7 keV band lags both lower energies (hard lags in the 
left panel) and higher energies (soft lags in the right panel). 

2. Different detectors show the same trends where the 6-7 keV 
band lags both lower and higher energy bands. 

This result is a generalization of the phenomenon shown in 
Figure 2 shows for a single flare in the. XMM-Newton light curve. 
It further confirms what has been seen in previously in shorter 
observations, where the iron line peaking at 6-7 keV lags the 
primary continuum that dominates energies at either side of the 
line. 
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Figure 4. Left: lag-energy spectra for MCG 5-23-16 at three frequency bins from a joint fitting of the XMM-Newton , Suzaku, and NuSTAR light curves. From left to 
right: 2 x 10 5 10 4 Hz, 10 4 2.5 x 1 0 Hz and 2.5 x 10 "-5 x 1 0 4 Hz. Right: the best-fit model for the three frequencies consisting of two Gaussian lines. 
Solid (blue), dotted (orange) and dashed (green) lines correspond to the three frequency bins in Figure 4 in increasing frequency, respectively. The normalization of 
the first Gaussian peaking at ~5 keV remain unchanged within statistical uncertainties, while the second Gaussian line peaking at ~6.5 keV decreased in strength as 
we move to higher frequencies. 

(A color version of this figure is available in the online journal.) 


To further study these lags, we present the energy-dependent 
lags. The previous study of MCG-5-23-16 showed that the lag- 
energy profile averaged over a relatively broad frequency band 
resembles a broad iron line. Our aim here is try to probe this 
profile for several frequency bands (i.e., different timescales). 

3.1.3. Detailed Energy-dependence of the Lag 

Light curves in 16 energy bins between 2 and 10 keV in steps 
of 0.5 keV were extracted. This was found to be an optimum 
trade-off between energy resolution and having reasonable 
uncertainties in the calculated lags. The lag of the light curves 
in each energy band was calculated taking the whole 2-10 keV 
band as a reference (excluding the band of interest each time 
so the noise remains uncorrelated). The frequency band where 
there is high signal was divided into three bins. The resulting 
lag-energy profiles are shown in Figure 4 (left). 

The lag-energy profiles are plotted for three frequency bands: 
2 x 10 — 5 — 10 — 4 7 10 — 4 — 2.5 x 1(T 4 and 2.5 x 10" 4 -5 x 1CT 4 Hz. 
For reference, the first XMM-Newton results in Zoghbi et al. 
(2013b) were for frequencies < 1 0 4 Hz, which corresponds to 
the first frequency bin. The results presented here and those 
in Zoghbi et al. (2013b) are remarkably similar despite the 
different observations, confirming that the observed delays 
are not transients, but rather a generic property of the object. The 
lowest frequency in the plots is 2 x 1 0 to a direct comparison 
with the Zoghbi et al. 2013b can be made. Including the lowest 
available frequencies makes the lags slightly larger (it can be 
seen in Figure 3) but the shape remain the same. 

The plots at the three frequency bins show clear line profiles. 
The shape of line appears to change with temporal frequency, 
peaking toward lower energies for the highest frequency bin 
(shortest timescales). The data above ~5 x 1 0 4 Hz are dom- 
inated by noise, and the lag-energy spectrum (not shown) is 
consistent with constant zero lag. The lag plot for all frequen- 
cies is the average of the three profiles, and it is dominated by 
the lowest frequency band. The energy profile in this case looks 
very similar to the lowest frequency profile. 

In all three spectra shown in Figure 4 (left), a constant or a 
smoothly increasing lag is ruled out with very high significance. 
Fitting these profiles with a simple model of a power-law and a 
Gaussian to characterize the existence and the general shape of 
the lag, we obtain the parameters shown in Table 1 (/ 2 = 49 for 
33 degrees of freedom). Although the Gaussian energies in the 
lowest two frequency bins differ only at the ~80% confidence 


Table 1 

One-Gaussian Model Fit Parameters 


Parameter 

Bin 1 

Bin 2 

Bin 3 

E (keV) 
cr (keV) 

7.06 ± 0.42 
2.13 ±0.42 

6.62 ± 0.45 
1.29 ±0.42 

5.03 ± 0.24 
0.49 ± 0.2 

Table 2 

Two-Gaussian Model Fit Parameters 

Parameter 

Bin 1 

Bin 2 

Bin 3 

E\ (keV) 
normi (s keV) 
E* (keV) 
norm 2 (s keV) 

4.9 ± 0.5 
174 ±78 
6.28 ±0.14 
1234 ± 280 

160 ± 59 
589 ± 223 

225 ± 80 
106 ± 120 


Note. * The energy was fixed between the three bins. 


level, the highest frequency bin is significantly (>99.999% 
confidence) shifted to lower energies. 

If instead we fit the spectra with a model consisting of 
two Gaussian lines, we find the best fitting model shown in 
Figure 4 (right). This model is motivated by the appearance of 
the lag profiles. In particular, the rise at ~4 keV appears to 
be persistent at all frequencies, while the lag value at 6-7 keV 
changes significantly. This model describes the data well with 
a x 1 = 39 for 35 degrees of freedom. The best fit values for the 
Gaussian lines are all consistent within ala uncertainty, except 
for the high energy Gaussian at the highest frequencies where 
it is not required. Therefore, we assumed the Gaussian lines to 
have the same energies and widths in all three frequency bins. 
This assumption allows the three bands to be directly compared. 

Although the model in Figure 4 (right) is simply a descrip- 
tive and not physical model, it provides valuable insights in 
understanding these lag-energy spectra. The low energy Gaus- 
sian line peaks at 4.9 ± 0.5 keV (90% confidence), and does not 
appear to change with temporal frequency (all normalization 
values are within ler uncertainty). The high energy Gaussian 
on the other hand, peaks at 6.28 ± 0.14 keV, and its strength 
decreases gradually toward higher temporal frequency (normal- 
ization: 1234 ± 280, 589 ± 223 and 106 ± 120 s keV for the 
three frequency bins in increasing frequency, respectively, see 
Table 2). 

If we further assume that the line profiles are due to a 
relativistic iron line, we can replace the two-Gaussian model 
with a relativistic line model. This does not take into account 
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Figure 5. Left: the best-fit model to the lag-energy spectra. The model consists 
of a power-law and a relativistic laor line. Solid (blue), dotted (orange), and 
dashed (green) lines correspond to the three frequency bins in Figure 4 in 
increasing frequency, respectively. Right: confidence levels for the outer radius 
of emission plotted as Ax 2 difference from the best fit / 2 value. 

(A color version of this figure is available in the online journal.) 



Energy (KeV) 

Figure 6. Peak of the light curve flare at different energies from the NuSTAR 
light curve. All the times are calculated relative to the first point. This is similar 
to the XMM-Newton example shown in Figure 2 (right). See Section 3.2.1. 


the full relativistic transfer function, but provides a tool toward 
understanding these lag-energy profiles. We use the laor (Laor 
1991) model, which is a delta function convoluted with a kernel 
of the Kerr metric. Its parameters in addition to the line energy 
and normalization are the inner and outer radii of emission 
(rin, r out ), emissivity index q and inclination. The model also fits 
the data reasonably well (x 2 /dof = 61/40). Individual spectra 
could not constrain all the model parameters separately, so we 
made the assumption that most parameters are the same for all 
frequency bins. r out was left as a free parameter and it is used as 
a proxy for the extent of the emission region. This will be the 
case when, to first order, the variability timescale (i.e., temporal 
frequency) is directly associated with the emission region size. 
The inclination of 38 deg and emissivity index of 2.2 from the 
spectral modeling in Zoghbi et al. (2013b) were used. 

The results of the fit are shown in Figure 5. The left panel 
shows the model for the three frequencies, while the right panel 
shows the confidence ranges for the outer radius parameter 
plotted as the A/ 2 difference from the minimum x 2 value. 
The figure clearly shows how the extent of the emission region 
inferred from the lag-energy profiles changes with temporal 
frequency, or equivalently with variability timescale. At the 
longest timescales (solid blue line in Figure 5), the emission 
region extends over a larger distance (>50r g ), and it becomes 
smaller when shorter timescales are probed (<5 r g for the 
shortest timescale). 

3.2. Compton Hump Reverberation 

After discussing the K-band reverberation (2-10 keV band), 
we present in this subsection the result of analyzing the NuSTAR 
light curves above 10 keV. Following similar steps to those 
used when discussing the iron K band, we start with the simple 
procedure then we consider a full frequency-resolved case. 

3.2.1. Simple Consideration 

The method in this section is similar to that in Section 3.1.1, 
where we fit a simple Gaussian line to an isolated flare in the light 
curve. We select the strongest flare in the NuSTAR light curve 
(~2.9 x 10 5 s in Figure 1) in twelve energy bands between 2 and 
80 keV. Light curves from the two modules FPMA and FPMB 
were combined. A plot of the flare peak time as a function 
of energy is plotted in Figure 6. All the points were shifted 
vertically with the same amount so that the first point is zero. 



Figure 7. Lag-energy spectra for MCG-5-23-16 using NuSTAR light curves. 
Left: lag-energy spectra for lags averaged over a wide frequency band covering 
6 x 1 0 6 -6 x 1 0 4 Hz. Right: lag-energy spectra at two frequency bands. Blue 
circles are for the same frequency band as the left panel. The orange squares 
are for frequencies 3 x 10 5 — 6 x 1 0 4 Hz. The central frequencies for the two 
bins are: 6 x 1 0 (blue circles) and 1 0 4 Hz (orange squares), respectively. 
(A color version of this figure is available in the online journal.) 

Although this procedure is again simplistic and does not take 
into account the stochastic nature of the light curves, it provides 
some useful insights. The peak of the flare seems to change 
with energy. Below 10 keV, the flare shifts are consistent with 
those in the XMM-Newton data and with the full frequency- 
resolved analysis of Section 3.1. The data above 10 keV show an 
increase with energy. Only one flare is used to produce Figure 6, 
and the features are statistically not very significant (a constant 
null hypothesis is not ruled out.). The following section studies 
the lag-energy profile systematically using full Fourier-resolved 
time lags. 

3.2.2. Fourier-resolved Time Lags 

The left panel of Figure 7 shows the lag-energy spectrum 
resulting from Fourier-resolved time lag analysis using the 
NuSTAR data. The light curves from the two modules FPMA and 
FPMB are fitted simultaneously. The left plot is for lags averaged 
over the whole frequency range covering 6 x 1 O ' 6 -6 x 1 0 4 Hz. 
The light curves were divided into segments of length ~30-40 ks 
with a 5 12 s time sampling. Eleven energy bins with logarithmic 
binning were chosen, 6 below 10 keV and 5 above it. 

Figure 7 (left) shows a general increase in the lag with energy 
with the ~6 keV feature that is consistent with what is seen 
in the combined data analysis of Section 3.1. There is also an 
increase above 10 keV with a possible peak at ~30 keV. The 
increase above 10 keV is very significant (e.g., against a constant 
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null hypothesis). The feature at ~6 keV is significant at >97% 
confidence when fitted with a Gaussian against a null hypothesis 
of a single power-law for the whole spectrum (this increases to 
more than 99.9% if the energy of the line is assumed known 
from photon spectral modeling for example). The turnover at 
40 keV is significant at the ~90% confidence when comparing 
a power-law fit with and without a cutoff energy. 

If our interpretation of the feature at 6 keV as due to relativistic 
reverberation in the iron K line (Section 3.1) is correct, then a 
reverberation in the Compton hump is expected above 10 keV. If 
the iron line is responding to the primary continuum variability, 
then so does the Compton hump produced as part of the 
reflection spectrum. The effect of a possible intrinsic continuum 
lag (i.e., lags in the primary X-ray source that has nothing to do 
with reflection) are discussed in Section 4. 

To investigate the frequency-dependence of the lag-energy 
spectrum in the NuSTAR data. Figure 7 (right) shows the same 
plot as in the left panel, and also a lag-energy spectrum at the 
highest frequencies (>3 x 1 0 5 Hz). This is slightly different 
from the plots in Figure 4. There, the signal was high enough 
to extract lag-energy spectra in three adjacent and independent 
frequency bins. Here, the two frequency bins overlap above 
3 x ICC 5 Hz, but they probe different ranges. A combination 
of large scatter and large uncertainties dominate the lag-energy 
spectra when independent frequency bins are considered. The 
central frequencies for the two bins considered in Figure 7 
are: 6 x 10~ 5 (blue circles) and 1 0 4 Hz (orange squares), 
respectively. 

Figure 7 (right) shows two key features. First, the lag-energy 
plot appears to be generally steeper at lower frequencies, due 
possibly to the contribution from intrinsic lags (see discussion 
in Section 4). Second, the two features at 6 and 40 keV appear to 
shift to lower energies at higher frequencies. The first statement 
is significant at more than the 4a level (for example we fit a linear 
model to the log of the energy we find slopes of 408 ±118 and 
75 ± 64 for the low and high frequencies, respectively), while 
the second statement is significant at ~2rr level (when fitted 
with two Gaussians, one for each peak, the energies shift from 
5.2 ± 0.5 and 42 ± 6 keV at low frequencies to 4.3 ± 0.7 and 
32 ± 4 keV at high frequencies for the low and high energy 
Gaussians, respectively). 

4. DISCUSSION 

This work presents an analysis that confirms the existence 
of reverberation lags in MCG-5-23-16. The lags are detected 
in new observations with XMM-Newton, Suzaku, and NuSTAR. 
Two further main results are presented: 1- The detection of a 
timescale-dependent change in the shape of the broad iron K line 
seen in the lag-energy spectra. 2- The detection of a time delay 
in the NuSTAR light curves that is consistent with reverberation 
in the Compton reflection hump. 

4.1. Reverberation below 10 keV: Disk Tomography 

The basic conclusion from the first result is that time lags 
seen in the short XMM-Newton data in Zoghbi et al. (2013b) are 
not only confirmed, but the new observations allow the emission 
from different regions to be separated. The results in the lag- 
energy spectra are model-independent as they are merely the 
outcome of comparing time shifts and delays between energy 
bands. The delays themselves have to be associated with a 
reflection process because there is a clear signature at the iron 
K energy at ~6.4 keV. 


The primary continuum as well as the reflection both con- 
tribute at all energies between 2 and 10 keV, but with different 
fractions at different energies depending on the shape of the re- 
flection spectrum (George & Fabian 1991; Ross & Fabian 1993). 
As discussed in Poutanen (2001; see also Zoghbi et al. 2013b), 
the lag-energy spectra measure the relative contribution of the 
reflection compared to the primary (power-law) component. In 
other words, the lag-energy profile measures the reflection frac- 
tion as a function of energy. For a simple spectral model con- 
sisting of a single power-law for the primary continuum and a 
reflected iron K line, as is the case here, the lag-energy spectrum 
roughly traces the shape of the iron K line itself. The measured 
shape for MCG-5-23-16 is clearly not a narrow line, indicating 
that relativistic effects are important. 

The analysis in Section 3.1.1 shows that by simply comparing 
flares at different energies, we find that the peak of a broad line 
lags energies below and above it. The widths of the flare changes 
accordingly. In a reverberation picture, the widths should 
be narrowest for the energies where the primary continuum 
dominates, and increase with increasing refection fraction. This 
is a simple expectation of a reverberating signal, where the 
reflected signal is a lagged, smeared , version of the irradiating 
signal, with the amount of smearing increasing with the size of 
the reverberating region. 

There is always a lag dilution factor caused by both primary 
and reflected components contributing to the energy bands of 
interest. The amount of dilution depends on the shape of the 
photon spectrum (e.g., Wilkins & Fabian 2013). Although the 
spectrum of the source in the observations used here will be 
discussed in detail in a separate work, the general shape is not 
very different from that reported in Zoghbi et al. (2013b) for 
the first XMM-Newton observation. The lags between bands 
measured here are of order hundreds to thousands seconds. 
The reflection fraction changes from ~30% at the peak of 
the relativistic line to ~10%-20% at the wings. This gives 
reflection-to-primary lags of order a few to several kilo-seconds. 
The light crossing time in seconds at a distance r (in units of 
r g = GM/c 2 ) from a black hole is GMr/c 3 = 50 Mir where 
M-i is the black hole mass in units of 10 7 M Q . The mass of 
MCG-5-23-16 is not known accurately, but it is of the order 
Mi ~ 1-5 (e.g., Ponti et al. 2012). Therefore, a delay of a few 
kiloseconds between the primary and the reflection, to within a 
factor accounting for the geometry, corresponds to a few to a 
few tens of gravitational radii. 

Interpreting the lag-energy spectra at different frequencies 
starts with the assumption that there is a correspondence 
between timescale and emission region size. Under this picture, 
different timescales (i.e., temporal frequencies) probe different 
emission region sizes. Combining this with the fact that lag- 
energy spectra measure the shape of the reflection component 
relative to the primary continuum that has a power-law shape, 
one can hope to measure the shape of the relativistic iron line 
emitted in different regions. The plot in Figure 4 shows this 
principle. The rest energy of the iron K line is 6.4 keV. The lag- 
energy spectrum at the lowest frequencies is broad and peaks 
at this energy, while the spectrum at higher frequencies peaks 
at lower energies. The line profile of a relativistic line from an 
accretion disk is well understood. Photons in the red wing of 
the line are emitted deep in the black hole potential, and they 
are emitted closer to the black hole. Photons emitted further 
from the black hole on the other hand, have a shallower black 
hole potential to escape and they are observed closer to the rest 
energy of the line. 
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An important point to note from the energy and frequency lag 
dependence is the fact that the red wing of the lag profile is the 
same at all timescales, while the core is more apparent at the 
longest timescales only, as expected from Fe K reverberation. 
As described in Cackett et al. (2014), the higher frequencies 
select against the longest lags producing a cut-out region in the 
core of the line that increases from low to high frequencies. 

The lag magnitude at the rest energy of the line changes from 
~700 s at long timescales to 0-200 s at the highest frequencies 
(Figure 4). This factor of ~4 change in the lag corresponds to 
a factor of ~7 change in timescale (the difference in frequency 
between the first and last bin in Figure 4), which is consistent 
knowing that the light crossing time (which scales with the 
lag) is proportional to the radius r while the timescale (most 
likely thermal in this case) scales with r ,/2 (Frank et al. 2002). 
We attempted to do a full modeling of the lag-energy spectra 
at the three frequencies using the simple lamp-post model in 
Cackett et al. (2014). We found an excellent fit to the lowest 
two frequencies with a source height of h ~ I Or,, . The model 
was unable to produce the highest frequency shape indicating 
possibly that the assumed geometry is not correct, and the corona 
is likely extended vertically. The observation seems to indicate 
that the shortest timescale response has to be emitted only at 
small radii illuminated by possibly a different source than the 
longest timescales. We defer exploring models with different 
illuminating geometries to a future work. 

4.2. Reverberation in the Compton Hump 

The analysis in Section 3.2 presents the first measurement of 
a time delay above 10 keV in an AGN. This is only possible 
thanks to the unprecedented hard X-rays sensitivity provided by 
NuSTAR. The existence of the lags in the data is significant and 
model-independent. 

Black hole binaries are known to have inter-band time 
lags above 10 keV (mostly from RXTE data; Nowak et al. 
1999; Grinberg et al. 2014), where a lag-energy plot shows 
a smooth increase with energy. Although their origin is not 
clear, they could be Comptonization lags (Kazanas et al. 1997) 
or propagating fluctuations lags (Kotov et al. 2001; Arevalo 
& Uttley 2006), with the latter being most likely. The key 
property in the Galactic black holes lags is their smoothness with 
energy. The absence of any features in their lag-energy spectra 
linked them to the primary Comptonizing continuum. The 
power-law shape of this Comptonization component explains 
the smoothness of the lag-energy spectra in Galactic black holes 
(Kotov et al. 2001). The question is: are the lags presented here 
for MCG-5-23-16 similar to those in Galactic black holes? 

Two points are considered in addressing the question. First, 
are the timescales comparable assuming black hole systems in 
AGNs are scaled-up version of Galactic black holes? Second: 
Do the lags have similar energy dependence? 

The lags in black hole binaries are measured at frequencies 
of ~1 Hz for Cygnus-Xl as an example (0.1-10 Hz; Nowak 
et al. 1999; Kotov et al. 2001). The ratio of the masses 
between Cygnus-Xl and MCG-5-23-16 is ~ 10/2 x 10 7 = 
5 x 10- 7 . Assuming the variability timescale is proportional 
to the mass (both viscous and thermal timescales scale with 
mass), a variability frequency of 1 Hz in Cygnus-Xl corresponds 
to a frequency of 5 x 1 0 7 Hz in MCG-5-23-16, which 
is smaller (i.e., timescale is longer) than probed with the 
current observations. Unlike the Galactic source timescales, the 
MCG-5-23-16 observations probe timescales that corresponds 


to thermal timescales a few and at most tens of gravitational 
radii from the black hole, hinting at a different origin. 

The second important point is the energy dependence. Lag- 
energy spectra in Galactic black hole are very smooth. Kotov 
et al. (2001) searched for any features in the lag-energy spectrum 
of Cygnus-Xl and ruled them out. The case of MCG-5-23-16 
is different. It is not just the NuSTAR data that shows a feature 
around the 6 keV at the same frequencies, but more importantly, 
the much better quality combined data from XMM-Newton and 
Suzaku presented in Section 3.1 do also. It is by looking at the 
same-frequency lag-energy spectra below 10 keV that the lags 
above 10 keV can be properly interpreted. 

Data below 10 keV strongly suggest a lag origin in relativistic 
reverberation. This comes from the existence of lags in the 
first place and then their energy dependence that resembles 
a relativistically broadened iron line. In this picture, energies 
above 10 keV are also expected to produce a lag feature 
corresponding to the Compton hump from the same reflection 
spectrum as the iron line. As in the case of the iron line, the 
increase in the lag-energy spectrum above 10 keV is produced 
by an increase in the reflection contribution in the spectrum. 
Furthermore, the analysis presented in Section 3.2 seems to 
indicate that there is a down-turn above 30 keV as expected 
after the Compton hump peaks and the reflection fraction starts 
to decrease again. Statistically, the turn-down is present at the 
2 a level when considering the whole frequency band. The 
significance is higher for higher frequencies suggesting that 
the turn-down is real. The contribution of the intrinsic continuum 
lags remains however unclear. 

Therefore, although continuum lags cannot be ruled out 
directly, the non-matching frequency scaling between Cygnus- 
Xl and MCG-5-23-16, the non-smooth lag-energy spectra in 
MCG-5-23-16 compared to Cygnus-Xl and the possible down- 
turn at ~30 keV, together indicate that the lags above 10 keV in 
MCG-5-23-16 are different from those in Galactic black holes. 
Lags in Galactic black holes are continuum lags characterized by 
a smooth increase with energy. They are due to Comptonization 
or more likely due to propagating fluctuations. The lags in 
MCG-5-23-16 are different, and there are strong indications 
that they are due to relativistic reverberation. 

The lag-energy spectrum of MCG-5-23- 1 6 using the NuSTAR 
data at two frequency bins (Figure 7) also indicates a generally 
steeper spectrum at lower frequencies. As we have indicated, 
continuum lags inferred from a simple mass-scaling of those 
in Cygnus-Xl are expected at 10 _7 -1CU 6 Hz, so there is a 
possibility they are contributing to the lowest frequencies in the 
bin giving a steeper profile. This means that lag-energy spectrum 
in this case is measuring the reflection fraction superimposed on 
top of a continuum lag that smoothly increase with energy. When 
the lowest frequencies containing the continuum are excluded 
from the frequency bin (orange squares in Figure 7 (right)), the 
general shape is now flatter, because it is only measuring the 
reverberation lags and the shape is essentially the reflection 
fraction only. This difference between low and high frequency 
energy-dependence is similar to what is seen below 10 keV in 
many objects (Zoghbi et al. 2010; Zoghbi & Fabian 2011). The 
case of Akn 564 for example (Kara et al. 2013b) shows clearly 
that the low frequency lags increase smoothly with energy, while 
at high frequencies relativistic reverberation traces the shape 
of a reflection component. With the Compton hump results in 
MCG-5-23-16, this property of two lag processes dominating 
at different frequencies is further supported by including hard 
X-rays. 
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5. SUMMARY 

This work presents three main results: 

1. Confirmation of iron K reverberation seen earlier MCG- 
5-23-16 using XMM-Newton data. New observations with 
XMM-Newton, Suzaku, and NuSTAR confirm the existence 
of the lags. 

2. Measurement the shape of the relativistic iron line at three 
timescales through the calculation of lag-energy spectra. 
The shape of the line changes with timescale. At long 
timescales, both the blue and red wings of the line are 
seen. At short timescales, only the red wing is seen. 

3. Using data from NuSTAR , we report the discovery of a time 
delay between energies >10 keV and the continuum. These 
lags are most likely due to reverberation in the reflection 
Compton hump. 

This work has been partly supported by NASA grant 
NNX14AF89G. The research in this article has made use of 
data obtained from the Suzaku satellite, a collaborative mis- 
sion between the space agencies of Japan (JAXA) and the USA 
(NASA). We made use of data from the NuSTAR mission, a 
project led by the California Institute of Technology, managed 
by the Jet Propulsion Laboratory, and funded by the National 
Aeronautics and Space Administration. Part of this work is based 
on observations obtained with XMM-Newton , an ESA science 
mission with instruments and contributions directly funded by 
ESA Member States and NASA. 
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